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The effective action for the charge density and the photon field is proposed as a generalization 
of the density functional. A simple definition is given for the density functional, as the functional 
Legendre transform of the generator functional of connected Green functions for the density and 
' the photon field, offering systematic approximation schemes. The leading order of the perturbation 

expansion reproduces the Hartree-Fock equation. A renormalization group motivated method is 
introduced to turn on the Coulomb interaction gradually and to find corrections to the Hartree- 
Fock and the Kohn-Sham schemes. 
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. I. INTRODUCTION 

CZ5 . 

The functional methods, based on the generator functionals for the different Green functions have proven to be 
useful in Statistical and Quantum Physics Q]. There are two kinds of important generator functionals characterizing 
the dynamics, corresponding to the connected and the irreducible Green functions. They are related by functional 
i-^j \ Legendre transformation and the latter is called the effective action 

Another widely used application of functionals is the density functional formalism In this method one intro- 

duces an external potential and studies the response of the charge density in the ground state Q. The key object of this 
, formalism is the density functional E gT [p], the ground state energy of the system with charge density distribution p(x). 

Few simple properties based on the Hohenberg-Kohn (HK) theorems [Q and the use of the Kohn-Sham (KS) scheme 
|^| obtained by means of the perturbation expansion, have led to powerful and widely applicable approximations. 
The problem left open is the determination of the density functional beyond the perturbation expansion. 

There are two further developments in the search for the appropriate density functional we relate in this paper. 
One is the so called adiabatic approximation to the density functional |sj combined with the coordinate scaling ||. 
In this scheme the Coulomb interaction is turned on gradually ]To[ | and the external potential is adjusted in such a 
manner that the ground state density remains invariant. The result is an interpolation, established between a non- 
interactive and the physical system. Another development is the realization that the HK theorems emerge naturally 
and trivially within the framework of the effective action formalism flTTJ . This relation may be found by means of 
the time-dependent variational method |l2| , or by considering the Helmholtz free energy for classical fluids p3| , as 
well. The perturbative emergence of the KS scheme |l4| is particularly simple in this language when the path integral 
formalism is used [ O . 

In defining the effective action we go beyond the scheme used in Refs. and jl4| by introducing an external 
source coupled to the photon field in order to simplify some intermediate stage of the computation. We generalize 
the adiabatic approximation within the framework of the effective action formalism by borrowing ideas from the 
q , functional renormalization group method. The renormalization group strategy was originally constructed to find the 
scale dependence in the UV region of QED and to treat the formal infinities in renormalizable Quantum Field 
Theories . It was realised later that the same method can be used to explain the singular properties of critical 
systems Further important applications include the partial resummation of the perturbation expansion Um, a 
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systematic setting of the operator product expansion jl9| and the solution of differential equations pOpl] . According 
to the original suggestion, one eliminates degrees of freedom in the Kadanoff- Wilson blocking and traces the evolution 
of the reduced dynamics as the system becomes simpler. A further generalization of the algorithm, the renormalization 
group method in the internal space allows us to turn on the fluctuations in the path integral gradually and follow 
the evolution of the dynamics in the meantime |2^ , |23|| . Such an implementation of the renormalization group idea is 
outlined here to perform a resummation of the perturbation expansion. 

The internal space renormalization group is the functional realization of the Hellmann-Feynman theorem |pi] , p5f for 
the ground state and makes possible to follow the change induced by a control parameter in the effective action. The 
resulting evolution equation is simpler than that of the adiabatic approximation. This is because the relation between 
the generator functional for the connected Green functions and the effective action, the Legendre transformation, 
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makes it possible to realize this scheme without the readjustment of the external potential, by turning on the Coulomb 
interaction alone. 

Another improvement offered by the use of the functional differential equation is the flexibility which allows us 
to turn the qualitative knowledge about the dynamics into a well chosen ansatz for the effective action. The local 
density approximation supplemented with the gradient expansion pfj is the usual ansatz for the density functional 
though excitations appear to bring in non-local contributions (^7). We propose in this paper another ansatz based 
on the idea of the multi-local structure of the density functional; it is similar to the cluster expansion in Statistical 
Physics. This ansatz avoids the limitations of the gradient expansion and gives a systematical truncation scheme for 
the functional differential equation. The solution of the evolution equation for the effective action which is truncated 
at the bi-local level, explored in the present paper, reproduces the HF scheme solution together with an infinite sum 
of radiative corrections to the exchange interaction and the ground state energy. The efficiency of the HF method 
might be explained by the fact that beyond these improvements all other radiative corrections are generated by at 
least tri-local, i.e. three-body interactions. 

A trivial lesson of the path integration method is that the exchange interaction originates from the one-loop level 
'zero-point fluctuations' of the photon field. The importance of this remark becomes evident when the KS scheme is 
sought and one finds that the renormalization group motivated resummation automatically generates the exchange 
interactions independently of the choice of the quasi-particle wave functions. In particular, an invariance property 
of the evolution equation with the effective action taken in the bi-local approximation makes possible to trade the 
exchange terms for a local potential energy and leads to the systematical improvement of the KS scheme. 

The organization of the paper is the following. We consider a system of electrons submerged in an external potential 
and subject of the Coulomb interaction. The generator functional W for connected Green functions of the density and 
the photon field and its Legendre transform T, the effective action are introduced in Section II. As the simplest example, 
Section III contains a short discussion of these functionals for non- interacting electrons. The most natural way to 
obtain the effective action is perturbation expansion in the Coulomb interaction. To perform a partial resummation 
of this series we introduce in Section IV a renormalization group motivated method, the evolution equation as the 
Coulomb interaction is turned on by starting with non-interacting electrons propagating in the physical external 
potential. It is shown in Section V that the leading order contributions of the perturbation expansion generate the 
HF energy functional. The main part of this paper is the discussion of the evolution equation for various truncation 
levels for the density functional in Section VI. Here we recover the traditional HF equations, their systematically 
improved version and the generalized KS equations. Section VI. G is devoted to a brief comparison of the effective 
action method with the density functional theory. The conclusions are given in Section VII. The Appendices present 
some details which are useful in the main text, the discussion of the physical interpretation of the effective action, 
the continuum limit of the Grassmann-integral representation of the ground state energy for fcrmionic systems, and 
Legendre transformation as a functional differential equation. 

II. GENERATOR FUNCTIONALS 

We consider in this paper fermions which have n s degenerate one-particle spin states and the system of N% a t = Nn s 
such particles propagating in the presence of a static external potential V ext and interacting with the Coulomb 
potential. After presenting the partition function which defines our system we turn to the functionals characterizing 
its dynamics, the generator functional for connected Green functions and its Legendre transform. We show finally 
that when Gauss' law, i.e. the Coulomb interaction is properly established for the external sources then the effective 
action simplifies in what it includes a special combination of the external sources only. 

A. Partition function 

The partition function 

Z = Tre-P<- H -»\ (1) 
serves as a projection onto the ground state as (3 = 1/T — > oo. The Hamiltonian H = Hi + H2 is given by 

H > = - 2 E L^ l '-'^^&^ (2) 

s,s' = l"' x > x 1 1 
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where J x — J d 3 x and we introduce the field variable 

oo 

^ s ,x = c S;W ^ n>x , (3) 

11=1 

with J \?* jj^l/m^ = <5n.m- The coefficient c s . n is taken to be the destruction operator, c S) „ = a s , n , for a single-particle 
state of the wave function $ n l in the operator representation and becomes a Grassmann variable in the path integral. 
The partition function is given in the latter case as 

Z = J V[^]V[^\D[u]e- s[ ^^' u \ (4) 

with 



(•5) 



We use the notation x = (x, t) and J x = J d 3 x Jq dt with the Euclidean time t £ [0, 0\ and the normalization 
J V[u]exp — \ f x (\7u) 2 — 1 is assumed for the photon functional integral. An UV cut-off was introduced to regulate 
the path integral and the normal ordering of H2, V^ xt — > V£ xt + e 2 /A XjX was performed in arriving at the action 
(|J). Note that ti is absent in the exponent of the Euclidean path integral and its only role is to control the non- 
commutativity of the kinetic and the potential energies. The Grassmann integral measure covers the trajectories 
in the functional space spanned by the one-particle wave functions ip s ,x.,t — Y) n Cs,n,t^n,x satisfying anti-symmetric 
boundary conditions, c s>n (f + [3) = —c Sjn (t), 

s n 

The range of n extends over a complete set of functions and the chemical potential /i is chosen in such a manner that 
the ground state contains the desired number of particles. For simplicity, we assume that all states are normalizable. 
When de- localized states are needed the index n has continuous spectrum, as well. The time-component of the photon 
field is purely imaginary for imaginary time as shown by the factor i in the minimal coupling. The chemical potential, 
the homogeneous component of the temporal component of the real time photon field remains real. 

Our basic assumption is that the formally e-independent one-body interaction Hi in Eq. (|J) includes the singular 
part of the interactions and the two-body interactions build up smoothly. We do not require the availability of the 
perturbation expansion in H 2 instead we assume that the true ground state, corresponding to the physical value of e, 
can be achieved smoothly by increasing e from zero in infinitesimal steps without encountering singularities. 



B. Connected Green functions 



The source-dependent generator functional for the connected Green functions is defined as 

P[M]e™ sTrlog[G_1+le " I "i [v(,I+ = ' j)1 ' [v(u+ = ' T)1+: ''' ( " + = ' ) . (7) 



where the inverse propagator 



Cr\, =,i, , ( (5t -|-A x - M + y x ext ) ,<s) 



together with the notation / • g = J f x g x , $x,x' — ^x,x'<^t,t') ^ x . x ' = S(x. — x'), and S ti t' = S(t — t') was introduced. 
Notice that the spatial gradient is diagonal in time, (V) X ,{ )X ',{' = # x>x '<5t,t' V x . 



3 



It is instructive to compare the role the field a and the Grassmann-valued sources 77, if play when the action 
S[ij)\ip\ — > S[ip^,ip] — icr ■ ipty + if • ip + ip< ■ 77 is used for electrons in the absence of the Coulomb interactions, 
e = 0. The dependence of the path integral on the sources 77, if is easy to find by diagonalising H\. On the contrary, 
the dependence on a cannot be obtained in a closed form. The sources 77, if control one-particle processes: the 
introduction or removal of particles, without influencing their states, the eigenvectors of H\. The source a leaves 
the particle number unchanged but modifies the states they occupy, induces a polarization, and its effect cannot be 
obtained in a simple form in terms of the eigenvectors of H\ . 

The imaginary time electron propagator 

Gx.x' = (0|T [ip^l] |0) (9) 
where |0) denotes the ground state is usually written in the spectral representation 

G x , x < = - £ / _*"2w- e " <(t "°" = - E (10) 

n=l w lL ° " 

for — /3 < t — t' < 0, /3 — > 00, where 



•— a x + y x oxt J ^ n , x = ^ n ^„ )X) (11) 

and J = J du/2'K. The set of wave functions <p n and the spectrum E n correspond to non-interacting quasi-particles 
in the presence of the external potential. We shall need a parameterization of the quasi-particles dressed by the 
Coulomb interaction. This will be achieved by writing the field operator (||) in terms of the dressed basis ^> n specified 
later, 

N 



G X . X ' = - ]T * n ^, iX/ (0|a n e-( H -"H*-*') a t,|o) 

n,n' — l 

N 00 
n— 1 n.n' — l 

for t > t' where the matrix elements of the time evolution operator are parameterized as 

- (0|a n e-( H -^4,|0) = 5 n , n .e-^-^+g n , n ,{t). (13) 

The relation 

(0\a n (H - n)a\\0) =E n - i x- d t g n , n (0) (14) 
together with the identification of E ni n > N with the energy increase when a particle is added requires 

Kg*,*® = n J%[ (15) 

The part n < N of this equation holds due to aJjO) = 0. Notice that g nn >(t) = 0(t), when (3 — > 00, assuming that 
the ground state |0) contains always N particles, i.e. the chemical potential is readjusted whenever the basis ^ n is 
modified. Finally, the functional trace can be written as 

TrA = f Ax,t,x,t = J2 I dtA n,t,n,t (16) 

J XJ n=l J 

where A nt t, n ',t' = / XjX ; (f>n, x -Ax,t,K.',f<f>n',x'- 

It is well-known that photon vertex functions display singularities when p — > 0. This corresponds to a qualitatively 
different dependence on the constant and the slowly varying components of the external sources. Notice that a given 
Fourier mode of the external sources, cr PtUJ , or j p>ul shifts or couples to the photon field u PiUJ , respectively. Since 
the restoring force of this photon field component to the vanishing equilibrium value is 0(p 2 ) the response of the 
system for the slowly varying external sources will be strong. In order to follow the smooth building up of the 
Coulomb interactions, we shall consider inhomogeneous external sources and their polarization effects only, i.e. we 
set <7 p= o, w = j P =o,u — in agreement with the absence of integration over the homogeneous, p = lu = modes of the 
photon field in QED §1. 
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C. Effective action 



The quantity of central importance is the effective action, the functional Legendre transform of the generator 
functional for connected Green functions. The reason this functional proves to be useful for strongly coupled particles is 
that (i) its minimum reproduces the free energy at finite temperature and the ground state energy at zero temperature, 
and (ii) its value at extended configurations gives the free energy or energy corresponding to a given constraint, cf. 
Appendix [A|, and finally there are well developed methods for its computation. The property (i) enables us to compute 
the ground state energy in the low temperature limit. The search for excited states is facilitated by feature (ii). 

The effective action is defined as 



r [P, v ] = -W[(T,j] + ia ■ p + j ■ v, 
where the induced charge density and photon field are given by 

SW[a,j] 5W[a,j] 



Px 



The inversion of these relations results in 



Sia x 
5p x 



3x 



Sv x 



(17) 
(18) 

(19) 



showing that the effective action has an extremum, obviously minimum in the ground state, j = a = 0. 

The free energy for a given density configuration p x and external potential v x , TT[p, v], is given by the functional 
integral 



p|- u j e « a Trlog[G- 1 +jeii]+i(«+iCT)-A-(u+iff)-io--p+J-(M+^( 

where the sources are found by inverting Gauss' law, Eq. (|l8|), 



(20) 



le-Px 











f m 


A (u x - 




+ 3x 









t ,n s Trlog[G" 1 +ieu] + i(u+io-)-A-(«+io-)+i-(u+icr) 



A 



(u x )a,j + ~cr a 



' Jx, 



V x = {Ux)a,i +-&X- 

e 

A particularly useful relation connecting the second functional derivatives of these functionals is 



S 2 W[a l3 ] 
8i<7 x 5ia x i 
S 2 W[a,j] 



S 2 W[*,j] 
Sia^Sj^i 
S 2 W[<7,j] 



s 2 r[ P ,v] 
i5 2 r[ P ,-u] 

<5iv<5p x " 



8 2 T[p,v] 
Sp x r8v m „ 

s 2 r[p,v] 

Sv i Sv n 



1 
1 



(21) 



(22) 



which can be obtained by evaluating the trivial functional derivatives like 



3 X ,x" — 



Spx 



S SW[a,j] 



8p X " 8p X " 5ia x 
6 2 W[a,j] S 2 T[p,v] , 



5 2 W[a,j] 8ia x , , S 2 W[a,j] 5 Jx 



Siax'Siax 6p x 
6 2 W[a,j] 6 2 T[p,v] 
Sj x 'Sia x 5p x "5v x ' 



5j x iSia x Sp x " 



(23) 



Sia X 'Sia x Sp x "Sp x 

etc. according to the chain rule. 

Expression (0) for the generator functional for the connected Green functions shows that the role of the test- 
potential v x is played by the source ia. The effective action T[p, v] gives the free energy of the system in the presence 
of such external sources which generate the charge density and photon field p and v, respectively, in the ground 
state. The photon dynamics is not traced by the density functional T[p] and the Legendre transformation yields the 
straightforward result 

E HK [p]=\\raJ-T[p\. (24) 

The change of variable a — > p is invertible so long the Jacobian Sp/Sa is non-singular. The latter is the connected two- 
point function which is positive definite in a non-degenerate ground state. The density functional is thus well-defined 
in a sufficiently small vicinity of the ground state density, p « p* . The only way the density functional may cease 
to exist is the singularity of the Jacobian Sa/Sp. According to the identity (p2]), this amounts to having a diverging 
eigenvalue for the density-density correlation function, the indication of an instability in the ground state. 
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D. Gauge invariant sources 



The two external sources, a and j, introduced in the generator functional (f7j) are not independent physically. In 
fact, the electrically charged environment of the system can be represented either by its external charge distribution, 
p ex t, or by its external electric field, w cx t which are related, Aii cxt — ep cx t = 0. This is Gauss' law which arises from 
the redundancy, the gauge invariance of the external sources. Since a/e and j appear as external electric field or 
density, Ac is identical with ej, in other words one expects no dependence on Ac — ej in W[a, j}. In a more formal 
manner, the generator functional must be invariant under the transformation 

a -> a - ex, j -> j + A%. (25) 

Such an invariance arises from the fact that there are two independent variables in Eq. (|^) but only one non-trivial 
functional, Trlog[G _1 + ieu — icr]. The relation (|l|) can be written as 

6W tot [a,j] _ 1_ A 6W tot [a,j] 
da e Sj 

where W tot [(T,j] = W[a, j] + j ■ ■ j showing that iy tot [tT, j) actually depends on the combination J = j only. 

The choice x — a i e an d X — ~~ A -1 • j gives 



W[a,j] = W 



0,j + -Aa 

e 



—z a ■ Act + -a ■ j = W 

2e z e 



j-^-j, (27) 



and the relation W tot [J] = W[J, 0]. It is easy to find the physical meaning of W tot [J]. In fact, according to gauge 
invariance or Gauss' law the external source j is coupled not only to the photon field u but has its contribution, 
A -1 • j to the photon field itself. The gauge invariant partition function of the electrons and external fields is 

e W tot [J] = j p[ u ] e n 3 Trlog[G- 1 + i e«~i<T]+^(u+J-i:)-A(«+i-j) 

D[w]e nsTl ' los[G_1+ ' ieu]+ ^ ( " + - J) ' A(u+ = J) . (28) 



The reduction of the effective action to a functional with a single variable can be achieved by rewriting Eq. (|2l]) as 

ST[p,v] 



8v 

whose solution is 



iep-Av, (29) 



T[p,v]=T[p,0]+iep-v--v-Av. (30) 

The three terms on the right hand side are the free energy for a given density and photon field as the sum of a 
mechanical term, the contribution of the electron orbits T[p, 0], the interaction between the charges and the photons 
iep ■ v and the free energy of photons — v ■ Av/2. We can construct the effective action for p, the density functional, 
by minimizing T[p, v] in v. The minimum is at fgrt/c] = ieA^ 1 ■ p, where 

r[p] = r[p,o]-i-p-i- P . (31) 

This relation indicates that T[p] is the complete free energy, that is the sum of the contribution of the electron 
structure, T[p, 0], representing the mechanical, correlation energy and the Coulomb energy of charge density p. The 
relation 

5T[p] _ ST\p,v] 

proves that the density functional is actually the Legendre transform of W[a, 0], or what will be useful for the 
computations below, of the generator functional of the complete electron-source system, W tot [a] , 

-T[p] +ia -p = W[a,0] = W tot [a]. (33) 

The reason to keep both sources a and j and work with the functionals W[a, j] and its Legendre transform 

e 2 1 1 
T[p, v] = T[p] + — p ■ — ■ p + iep ■ v - -v ■ Av (34) 

instead of using a simpler single-variable functional will be explained in the next section. 



G 



III. NON-INTERACTING ELECTRONS 



There are two different problems intermingled in computing the density functional T[p\. One is the expression of 
the external source, the external potential of the density functional theory, a in terms of the ground state density, p 
when the Legendre transformation is appl ied to the non-interacting electron problem. We shall explore this problem 
in this section. We show in Appendix C 2 how to obtain the same result by solving a functional differential equation. 
The other problem, left for the subsequent section, is to take into account the Coulomb interaction. 

We start with the generator functional for non-interacting electrons 



W free [a] = nsTrlogQ^ 1 = ^TrlogG" 1 - n s V — Tr(G • ia) n (35) 

z — ' n 

n—l 

was introduced. Here a is considered as a diagonal operator acting on the 



where the notation Q x 1 



space-time <T Xj t, X ',t' = iJx.x'^.t'^x,*, (G ■ a) x ,x' — G x _ x ia x . It is easy to understand the role of the infinite series in this 
equation. The free energy of non-correlated particles in the presence of an (Euclidean) external electric potential ia/e 
is p ■ ia, where p is the density. This is the term n = 1 in Eq. (j35"[). We have to correct this expression if the particles 
are correlated at different space-time points. In order to take into account the influence say of a particle at point a 
on the probability of finding another particle at point b, we add a term ia ■ ■ ia to the free energy where Cf£\ 



is the connected density-density correlation function. The contribution with a given n in Eq. 
rt-particle correlations in the free energy. 

For the effective action we need the inversion of the equation 



corrects for the 



TisQ x , x 



n=0 



According to Feynman rules, the external source a is attached to graphs by the particle-hole propagator 



Gx,y — 



77>sG X ^yGy_ X . 



(36) 



(37) 



Therefore one expects that p, the Legendre-pair of ia, will appear in the graph in the combination r = G 1 • (p — p*). 
Since the charge density is non-vanishing in the ground state, 



we look for ia as a functional power series in p 
form 



P x 
-P* 



Tl s G xx 

In this manner the inversion of Eq. 



(38) 

|) will be sought in the 



E 

n=l Jai < 



1 



i free 



I ^ l a,ai,...,o„ ' ai 



(39) 



The index pairs (x a ,t a ) are denoted by Latin indices a, etc. By inserting this expansion and Eq. (35) for Q into Eq. 
(Bq) and comparing the terms with identical order in p on both sides up to C(r 3 ) we find 



— $a,b, ^q^c — 2n s / G a L d Sd,b,c, ^aXc,d ~ 3! 



1 free 



\ free 



G~ 



Se,b,c,d ~\~ 2fl s 



f,g 



with 



Sa 1 ,a 2 ,...,a rl — , ^ ] ^op(i) ,Op(2) ^"p(2) ,«p(3) ^"p(n) ,«p(l) 



(40) 



(41) 



Pes n 



where the sum extends over all permutations of n indices. These expressions for the free gas correspond to the zeroth 
order ones in e for the interacting gas. In the last step we insert ( p9[ ) in the definition of T, and rewrite the result in 
terms of the matrices 



: I da 1 ,a 2 ,...,a n '~ r a 1 , Xl lj 'a2.X2 ' ' ' ^Qn,; 



(42) 



as 
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oo 

r free [ P ] = -vK free [(r] -(- id • p = 4 ree + £ / -r£ ee . iar 

n=2 J«i,..,a„ n - 



{P~ P*)a x •■■{p-P*)a n 



(43) 



with 



fc- = -n s TrlogG-\ r^ e = G^ b , T%^ c = 2S aAc , T^ d = 3\S a , b<c , d + 2 • 3! / S^G^S/^. (44) 



The interpretation of the form (43) goes along the same line as those of the free energy (p5[), namely T^ 66 with 
k < n different space-time point arguments represents the contribution to the free energy due to the fc-particle 
correlations when the particle density is modified compared to the ground state. 
At the minimum p = p* we find the ground state energy 



E„ 



lim — n s TrlogG 1 , 

/3^oo p 



(45) 



by assuming that the matrix G is positive semidefinite. 



IV. EVOLUTION EQUATION 

We tackle the second problem mentioned above, of the Coulomb interaction by a method motivated by the renor- 
malization group. For this end we introduce a parameter A by modifying the action S — » S\ 



(46) 



For A << 1 the Coulomb interaction is negligible and the real physical problem is recovered at A = 1 because A 
controls the photon field fluctuations. After the rescaling u — > Xu we find the action 



(47) 



which can be obtained from Eq. (^|) by the replacement e — > Ae. Note that the source term contribution in the 
generator functional, (0) is not scaled for simplicity. 
The functional differential equation, 



d\T x [p, v] = T[T X \p, v];X,p,v] 



(48) 



obtained below for the effective action of the model with a given control parameter value will then be integrated to 
find the true effective action, 



T[p,v] =T 1 [p,v] =T\ [p,v] + / d\T[T x [p,v]; \,p,v] 

JX 



(49) 



For a sufficiently small initial value, Ao * 0, the initial system is perturbative and T\ can be replaced by the one-loop 
expression (f78|). 

The evolution equation (ktq) can be obtained by noting 



d x T x \p,v} = -—Wx[a,j] 



iee -w x [aj] J x>[^]v[^}V[u] J ^ x u x e~ Sx ^ tfl^+i-* 

J X 



(50) 
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where the subscript c in the expectation value denotes the connected part of the appropriate Green function. We 
express the connected Green function by means of the functional derivative of the generator functional W[a, j] and 
find 



d\T x [p, v] = ie 



6 2 W x [a,j] 6W x [a,j]6W x [<T,j] 



Eq. (|2^) is then used to arrive at a closed differential equation 



8ji 



(51) 



d\Tx[p,v] = ie 



Tr 



S 2 Tx[p,v] 
SpSv 



+ ip ■ v 



The reason to keep the dependence on the sources a and j is to arrive at the evolution equation ( p0| ) with the 
hand side that can be generated by the second functional derivative of the generator functional W. But Eq. 
allows us to reduce the problem to functionals of a single variable. The inverse of the operator 



<5 2 r A [ P ,*»] 

SpSp 

s 2 r x [ P ,v] 

SvSp 



s 2 r*[ P M 

SpSv 
SvSv 



a 2 r A [ P ,o] 

SpSp 

iXe 



iXe 
-A 



(52) 



r (fof 



(53) 



S 2 W x [a,j] 
SiaSia 

S 2 Wx[<r,j] 
SjSia 



S 2 W x [a,j] 

SiaSj 
S 2 W x [<r,j] 

Sjsj 



-A 
—iXe 



a 2 r A [ P ,o] 

SpSp 

a 2 r A [p,o] 

SpSp 
-1 



A + A 2 e 2 
A + A 2 e 2 



s 2 r x [ P ] 

SpSp 



iXe 



—iXe 
^ 2 r A [p,o] 

SpSp 



a 2 r A [ P ] 

SpSp 



A <5 2 r A [ P ,o] 



A 



iXei 



s 2 rM 

SpSp 



A 



A 2 e 2 



i_ ( s 2 r x [ P ] \ 

A ' V SpSp J 



SpSp 



A 



+ A 2 e- 



(54) 



according to Eqs. 



and (34) and the evolution equation can be written in the form 



dxTx[p] = -Ae 2 



Tr 



S 2 Tx[p] 
SpSp 



1 

A 



(55) 



in terms of the density functional Tx [p] ■ 
it. By writing the right hand side as 



Note that this is an exact equation, no small parameter was used in deriving 



dxTx[p] = 



x,y,t 



Ae 2 



4vr|x - y| 



= Aa 



(V , i,i'0x,tV'y,t'0y,t>' 



x,y,t 



x - y| 



(56) 



where a = e 2 /Air, we find that the evolution equation for T[p] follows the accumulation of the free energy of the 
electronic structure in the absence of external photon field. The evolution equation for T[p, v] is more complicated 
because it takes into account an arbitrary, externally prescribed photon field. 

There are two kinds of contributions to the evolution equation: First, the explicit A-dependence in ( f47j ) gives the 
resummation of the perturbative contributions of the interactions between the particles. Since we assumed that the 
e-dependence of the ground state is continuous we may have to adjust the chemical potential p, — > px in order to 
keep the particle number iV tot fixed in the ground state. The chemical potential dependence is step function like 
for [3 — > co and the contribution to the evolution equation containing dxpx is made vanishing by the adjustment. 
Another, implicit A-dependence is generated by the following optimization of the choice of quasi-particles. The density 
functional will be minimized with respect to the one-particle basis $„^ x to take into account the interaction with the 
external potential V ext . This variational step will be carried out at each value of A and yields another contribution 
to the evolution. This second A-dependence which influences mainly the ground state properties can be obtained by 
performing the perturbation expansion of the HF equations in AA and retaining the leading order terms. In this 
manner the usual iteration of the HF equation which is supposed to minimize the (one-loop) energy functional will 
be carried out paralel to the resummation of the loop expansion. The resulting scheme maps the initial conditions, 
the non- interacting one-particle wave functions into another set. This latter is used to write the exact resolvent ( p^| ) 
in the spectral representation. 



9 



V. PERTURB ATIVE TWO-BODY FORCES 



Our goal in this section is to obtain the effective action ([17]) up to the terms quadratic in the fields by first computing 
the generator functional (|t]) for the connected Green functions and after then performing the Legendre transformation 
( |i"8| ) . For this end we expand the fermion determinant of the third line of Eq. . 

We do not follow the perturbation expansion in e, instead first we integrate out, at least formally, the electron field 
and after then apply the saddle point expansion for the resulting path integral over the photon field. On the one 
hand, the improvement in following this strategy is that one can retain an infinite subset of Feynman graphs arising 
from the perturbation expansion in e. On the other hand, the non-systematic nature of this scheme is that there is no 
small parameter for the saddle point approximation of the photon path integral. The remedy of this latter problem 
is to apply a cluster-expansion like approximation scheme for the photon path integral, to truncate it by ignoring the 
n-body correlations with a certain n. 

First we recover the HF energy functional in the leading order, 0(e 2 ) when the two-body correlations are retained. 
After that we point out that the direct Coulomb and the exchange interactions correspond to the leading order terms 
in the scheme where the two- and three-body correlations are ignored, respectively. Another point hinting at the 
different dynamical origin of the direct and the exchange terms is that they correspond to the tree and the one-loop 
level evaluation of the photon path integral which indicates that they describe correlations due to strong (classical) 
electric field and due to correlated quantum fluctuations, respectively. 



A. Tree-level contribution 

Let us rewrite the generator functional for connected Green functions given by Eq. (0) as 

e ff|,j| = J V[u]e- S ^ l[u] (57) 
by expanding the logarithmic function in the action S c ff [u] of the effective photon theory, 
S eS [u] = -W lree [-eu] - i (u+ V) • A • (u + ^-aj -j ■ (u+ \r 

= -nsTrlogG- 1 + n s Tr ^ ^— ^- (G ■ ieuf - ^ (u + ^ • A • (u + ~tr^ - j ■ (u + -t^ . (58) 

We consider Eq. ( |57| ) as the definition of an effective theory for the photon field in the presence of the electron system. 
Similarly to the interpretation of (|35"1), the infinite series in the effective action corrects for the correlations between 
the electrons when propagating on a given external background. 

The loop expansions of the original functional integral (0) and of the effective theory (|57]) are not equivalent because 
the trace in the exponent of ( p7| ) comes from one- loop order of the Grassmann integral in (^|) . Nevertheless one can 
apply the saddle point expansion to the effective theory. On the tree-level of this scheme we approximate the path 
integral by the maximum of the integrand, W[a, j] w W tTCC [er, j] , where VF tree [a, j] is the value of the effective action 
taken at its minimum u = u c \. The saddle point u c \ satisfies the equation 



= „ iens g ( _i r [(G . ieUclT G]xx _ 



/;=() 



A / 1 S 

A • (u c i + -a) 

e 



- ix = 0. (59) 



When the terms 0((eu c \) 3 ) in the exponent of (J57]) beyond the polarization of the ground state are neglected then we 
find 

ltd = D ■ (j - iep* + ^A • a^j , (60) 

where 

D~ l = -A + e 2 G (61) 

stands for the photon propagator. We can write the corresponding electric energy by introducing the induced electric 
field in the ground state, E = Vw c i, as 
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2 Su \ 5uSu 

= n s TrlogG _1 + -Ar<7 ■ A ■ a +-a ■ j + - [ i - iep* + -a ■ A ) • D ■ ( j - iep* + -A ■ a 
2e 2 e 2 \ e / \ e 

= n.TrlogG- 1 + i J E 2 + ^ Mcl . G • u d + -^a ■ A ■ a + X -a ■ j + G(e 4 ). (62) 

Here the first term in the last line is the non-interacting contribution to the mechanical energy of the electron system, 
the second and third terms represent the energy of the photon field u c i and of the polarization, respectively, whereas 
the last two terms account for the energy of the external sources. 



B. HF energy 



The next step is to consider the path integral ( p7| ) in the one-loop approximation and retaining the two-particle 
correlations only. First we take the case a = 0, when we find 

1 



W[0,j) = W ucc [0,j] + W 1 -^ = |(i - iep*) ■ D ■ (j - iep*) + W (63) 
where W = nsTrlogG" 1 + VF 1- ' 00 * 5 and the one-loop correction W 1 loop is independent of the source j, 

W^ loop = -^Trlog ^f c f [ ° ] = -^TrlGgi?- 1 + ~Trlog-A. (64) 
2 ouou 2 2 

We shall use the parameterization 



Y? h = -n e 2 



N 



N 



,n,m— 1 



■EEv(t-^ E " (H ' ) ^,^;,A,xt;y 

m— 1 n,n' 
N 

E E sm,r,'(<-<>- J5 » (t - t ' ) * m ,^;- J[ f„,x$; ![) 

n— 1 m,m' 

E 5n,n'(*-*')fl r m ) m'(*-*')*T.,x*^ ) x'*m 1 x'*^/ 
n.n'm.m' 



(65) 



for the photon self-energy due to the particle-hole polarizations YP h = —n s e 2 G for —j3 < t — t' < 0. The complete 
generator functional is 



W[<r,j] = -^(T • (G - e 2 G • D ■ G) ■ a + iea ■ p* + -J ■ a + ^{j - iep*) ■ D ■ (j - iep*) + W 



(66) 



according to Eq. ( J27 

The leading order, 0{e 2 ) HF energy of the ground state corresponds to<7 = - 7 = 0as/3— > oo, 



where the last term comes from 



I 1 °° 1 
Trlog-D- 1 - - = -E 



A 2 " n 

n=l 



; 2 G- 



(67) 



(68) 



The first two terms in Eq. (]67]) stand for the expectation value of the one-body Hamiltonian Hi in the state where 
n s particles are found with wave function ^ nx , n = 1, • • • ,N. The third term represents the exchange energy, 
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Thus, one obtains 



„, i J X 

TV 

E 



n=l ' 

n s e 2 



h 2 

-— A + F x ext 
2m 



1 



2 2 N 

nze 



E 



n, m— 1 x,x 



'47T X - X' 



47T X — X 



(70) 



The different powers of n s in front of the direct and the exchange terms remind us of their different origins. The 
multiplicative factors 0(n 2 ) and 0(n s ) characterize the contributions of the fermionic and bosonic determinants, 
respectively. 



C. Density functional 



We re-derive the HF equations in this section by following the Legendre transformation ( |l7| ) and minimizing the 
resulting energy density. According to (§7)), the density functional satisfies the equation 



- T[p] + ia ■ p = W 



0,-Aa 



(71) 



Using the one- loop result, Eq. 



the source a is to be eliminated in favor of the density p by means of the equation 



ip = A-D- 



1 



r A(T - ip" 



1 



Aa + 0(e 2 ) + 0{a 2 ). 



(72) 



One has to be careful in determining the orders which are reliable because of the factor 1/e appearing on the right 
hand side of Eq. (|7l|) in the argument of W . Due to this factor the higher order contributions in a will not be 
suppressed by e. As a result we have a double expansion in e and a. By inverting these relations one should use the 
perturbative expression 



£ = -l- e 2 i.G-i + 0( e 4 ) 
A A A 

for the photon propagator as late in the computation as possible. Instead, by relying on ( |6l| ) we find 

ia = G- 1 ■ (p - p*) + 0(e 2 ) + 0{p 2 ) 

and 

m = \p-G- 1 -p+ \p* .Q-i.f-p. g- l p* -W + 0{a 2 ) + G(e 4 ) 



(73) 



(74) 



(75) 



(76) 



(77) 



showing that the density fluctuations around p* generate two kinds of energy increases in the electron system: the 
contributions corresponding to the terms G _1 and e 2 /A are related to the mechanical correlations and Coulomb 
interactions, respectively. 

The complete effective action is found from Eq. (34), 



e 2 ~ e 2 

Q = G-G- — -G + 0(e 4 



where the one-loop electron-hole inverse propagator Q 1 has been introduced, 

s- = «}-'. A , 

It is instructive to compute the effective action without the photon contributions, 

r[p, o] = \{ P P *) ■ (g- 1 + £ )■ (p - p*) - W Q + 0(a 2 ) + G(e 4 ), 
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T[p,v] = -n.TrlogG- 1 + i(p - p*) ■ g- 1 ■ (p - p*) + iep ■ v - ]-v ■ Av + ^Li-logi?- 1 - ^Trlog-A. (78) 

The first term is the sum of the single particle energies at p = p* , the second term represents the change of the 
polarization energy due to the excitation of the electron system, the third one describes the interaction of the charge 
density with the electric field, the fourth one is the energy of the electric field, and before the last term contains the 
exchange energy. 

The ground state energy functional is found in the low temperature limit as 



1. 

/3™<i) (3' 



E gI [p,v] = Jim -T[p,v}. (79) 



The minimum of the density functional is reached at 



Pgr = ( 1 + e 2 G- j^-j ■ p\ v gT = iej^- ■ p gv . (80) 



Inserting these in (|79|), one obtains the ground state energy ( [70[ ) in the leading order. Hereby, the first term on the 
r.h.s. of (|78|), the fcrmion determinant provides the sum of the single fermion energies, the second term does not 
contribute at 0(e 2 ), the third and fourth terms give rise the Coulomb-energy, whereas the exchange energy (fifth 
term) occurs due to the boson determinant. The charge density of the ground state is modified with respect to the 
value p*, but this does not influence the energy of the ground state in the order 0(e 2 ). 

The second functional derivatives of the effective action give the inverse propagators. One might wonder why do 
the leading order inverse propagators appear in the one-loop level effective action ([78|). The inverse of the free photon 
propagator, —A, for example does not contain the polarization of the electron ground state. But the latter can be 
taken into account by minimizing in the density for a given external photon field configuration. In fact, for a given v 
the minimum of ( |78| ) is reached at p gr [v] — ~ieG ■ v + p* and the effective action for the photon field only, 

T[p gI [v},v] = -?i s TrlogG" 1 + -v ■ (-A + e 2 G) ■ v + iep* ■ v + ^Tr log D' 1 - - Trlog-A (81) 

reproduces the one-loop inverse propagator. This is the reason of the absence of the photon self-energy in 1/A on the 
right hand side of the evolution equation (|55|). Similar consideration gives T[p, v gI [p\] = T[p]. 



D. Exchange-correlation energy 



It is worthwhile noting the elusive source of the exchange interaction. In the usual, variational derivation of the HF 
equations the direct Coulomb and the exchange terms correspond to the diagonal and the cross terms in computing the 
expectation value of the Coulomb potential energy operator sandwiched between Slater determinants. In the second 
quantized formalism the exchange term comes neither from the fermion determinant which incorporates the single- 
electron aspects of the dynamics nor from the Slater determinant because the states in the Fock-space are properly 
(anti)symmetrized by construction, without using any determinant based on non-symmetrical wave functions. Instead 
it arises from the photon functional determinant, showing its genuine quantum nature. 

According to the remarks made after Eq. (|||), the fluctuating photon field u induces free ener gy c hanges 0(u n ) 
through the n-particle density correlations, the term corresponding to n of the exponent in Eq. (p7[). On the tree 
level ( p2[ ) we find the interplay of the 0(u), local and the 0{u 2 ) two-electron correlation contributions. They make 
up the direct Coulomb contribution of the external sources and their polarization effects. The common feature of 
graphs appearing in the tree-level solution of the effective theory is that the external legs end in the combination 
j — iep* — Aer/e, the inhomogeneous term in the equation of motion, (|59|). The external legs are attached to the 
rest of the graph by the photon propagator, D. These contributions describe the Coulomb energy of the local charge 
distribution and the polarization induced by the external local classical sources, a and j and internal local fluctuations, 
P*. _ 

The one- loop contribution, W 1 loop , is the 'zero point fluctuation' of the photon field. The normal modes of the 
photon field in ( |57j ) are coupled to the 0(u 2 ) two-particle correlators and VF 1-10013 is the sum of these contributions 
to the energy. The n-th order contribution of this correlation in the perturbation expansion corresponds to the term 
n in the expansion (|6^) where according to Wick's theorem the In insertions of the photon field are broken up into n 
photon propagators connecting the density correlation functions. The common element of the one-loop level graphs 
is that they have no external legs since a, j and p* couple to the first power of the photon field (minimal coupling) 
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and they drop out from the second functional derivative of the effective action. These contributions represent the 
Coulomb energy of the bi-local or higher order correlated fluctuations, the leading order contribution being the bi-local 
'exchange densities' introduced in Ref. [ ^0| . 

We considered above local densities and two-particle correlations. When higher order correlations are retained in 
the equation of motion ( |59| ) to determine the saddle point then higher order density correlations appear on the tree 
and the one-loop level, as well. 

The traditional direct Coulomb and exchange terms are the leading order C(e 2 ) contributions of the tree and the 
one-loop level graphs, respectively. These and other, higher order graphs of the tree and the one-loop sectors can be 
related to each other by 'opening' the electron lines of p* at the leg of a tree-level graph and attaching to another 
electron line, obtained in a similar manner. This explains that despite their difference in the effects they incorporate 
these graphs happen to be in the same order of e and h (when restored for real time). This may happen because the 
effective action (57) contains the sum of the tree and one-loop level contributions of the Grassmann integration and 
there is no further small parameter for the saddle point expansion within the effective theory. 

The exchange, maybe better to say correlation interaction nevertheless makes its appearance even on the tree-level 
of the effective theory, i.e. in the classical approximation for the photon field. In fact, the term 0(e 2 ) in the second 
equation of (p3) can be interpreted as the contribution of the bi-local exchange density to the electric energy. 



VI. MULTI-LOCAL APPROXIMATIONS 



We shall use in this section the evolution equation to perform a partial resummation of the perturbation expansion 
in the Coulomb interaction. The non-perturbative external potential is taken into account variationally by minimizing 
the effective action against the choice of the one-particle wave functions $f n ,x at each stage of the evolution. The 
evolution equation will be studied in different truncation schemes. For systems with homogeneous ground state the 
expansion in the inhomogeneity of the excitation, the gradient expansion, is a generally employed truncation scheme of 
the effective action. The localized bound states render this method useless and we have to turn to another classification 
scheme, alluded to in the previous section, the cluster expansion. 

Each truncation scheme S will be characterized by a subspace As of the space A of effective action functionals and 
the evolution equation is approximated by projecting it into As- In order to make this projection easy to define, the 
subspace As must be linear. Terms like J x ^ x Jxi,---,x e Pxl ' ' ' Px e e m the effective action represent the interactions 
between I clusters and will be called ^-local. A functional will be called £-local in general if it contains €-local clusters 
only. The linear subspace made up by the sum of n- local functionals with n < £ will be denoted by Cg. The effective 
action, T[p], contains multi-local terms with arbitrary high values of £. The approximation fi where terms up to l- 
locality are retained in a free, unconstrained manner corresponds to projecting the effective action onto the subspace 
Af t =Ct- 

It will be important to preserve certain functional structures within a given locality subspace. For this end we refine 
our scheme and introduce the constrained £-local approximation, a, based on the subspace A Ce of £-local functionals. 
Some complications arise form the fact that the functionals obtained by imposing a particular parameterization 
generally do not span a linear space and the restriction of the evolution equation for such a space is not uniquely 
defined. 

The linear subspace of the mixed scheme fkCg with k < t will be constructed by considering the functionals which 
are the sum of two functionals. The first one is a constrained one (q). The other belongs to the unconstrained class 
(fk) but chosen in such a manner that it does not interfere with the previous component. This will be achieved 
by orthogonalizing the two components. In such a manner the subspace of the mixed scheme is defined as Af kCe = 
A Ce + (A f k n Aj: ) where Ajr e denotes the orthogonal complement of A Cl . We want to avoid the detailed definition 
of the functional spaces in question, it is sufficient to say at this point that Cg is spanned by functionals of the form 
Jxi x "fx 1 ,---,xkPxl ' ' ' Pxt an d the scalar product is chosen in such a manner that the functionals p" and p™ be 
orthogonal for n ^ m. 

Once the subspace of the approximation has been identified, both sides of the evolution equation ( ]48| ) is projected 
onto Af e or Af k r\A^r in the scheme fg or fkCg, respectively. The evolution of the constrained component in the latter 
case is determined by imposing the evolution equation at certain suitably chosen configurations p x . 

One may distinguish three levels of the description of the ground state of N particles. 

1. The ground state is characterised by a wave function of 3N variables. 

2. The spectral representation is evoked to parametrize the Green functions, it is based on quasi-particle wave 
functions, N of them corresponding to the filled states and the rest for the virtual excitations. 
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3. A single function, the density is used and the coefficient functions of the constrained scheme are treated on a 
space-time lattice without any reference to wave functions. 

The first level is for the first quantized formalism, the second exploites the simplifications inherent in second quantiza- 
tion. These levels are certainly justified for small systems but lead to complications when the ground state of a large 
number of particles is sought. In fact, the computational time is at least 0(N 4 ) for a C2 scheme at description level 2 
(Hartree-Fock) and definitely worse for level 1. Level 3 offers more effective algorithms for large systems. Namely, after 
the diagonalization of a non-interacting Hamitonian for the initial condition, a 0(N 3 ) step, the Coulomb interaction 
can be taken into account in an f\C2 scheme. This latter is based on the density Green function, a space-time matrix 
of dimension N p . The power < p < 1 depends on the resolution in the space-time needed to reach the desired 
accuracy. In fact, for a fixed, particle number independent resolution we have p = and one finds p = 1 when the 
dimension of the one-particle sector of the Fock space is kept proportional to the particle number. The computational 
time is therefore 0{N 2p ), slower increasing with N than for the diagonalisation of the non-interacting Hamiltonian 
in the initial condition. 



An /1C2 scheme of level 3 is presented first in Section VI B , followed by a number of approximations of level 2. A 
simple C2 scheme where we neglect the terms C(e 4 ) and restrict ourselves to the parametrized bi-local approximation 
reproduces the HF equations. A possible simple improvement is an /1C2 scheme which may go beyond the leading 
order of the perturbation expansion. The most involved scheme sketched in this paper is f2- Finally we mention a 
modified f\Ci approximation leading to a generalization of the KS scheme. 



A. Quasi- particles 

The internal structure of an extended state is characterized in perturbative Quantum Field Theory by form factors. 
In the functional setting quadratic functionals usually correspond to free point-like particles and the form factors 
are hidden in the non-quadratic, higher order terms. This situation can be changed by considering the propagator, 
the kernel of the quadratic part as a free variable jy], occuring on the description level 3. We choose here a simpler 
short-cut for the level 2 description by parameterizing the quadratic kernel by means of one-particle wave functions, 
as given in Eq. ([To]) . This expression will serve as the definition of the quasi-particles making up the ground state 
and the excitations above it. 

The success of the effective action as a quadratic functional in the density suggests the form 

r[ P ] = \{p- P *)- Q- 1 ■ (p - p*) - ^f P ■ ^ ■ p + c[ P ] (82) 

where the 'potential energy' C[p] is introduced to take into account the higher order radiative corrections, Q and p* 
are given by (|7^) and (^), respectively. We write the evolution equation imposed for a fixed, A-independent p as 



\(P P*) • d x g~ l ■ (p - p*) - d xP * ■ g- 1 ■ (p - p*) + d x C[p] = -Ae 2 £(-l)"Tr Q ■ (c^ ■ 0) 



2 

where the notation 



n=0 



x\ n 1 



(83) 



Kl+C?l (84) 



S 2 T 

SpaSpb 

is introduced. Note that the summation on the right hand side of (^3|) can be constrained to n < I when the effective 
action is sought in £-local form because contains ultra-local terms in the density. The evolution equation ( |83| ) will 
be considered in two different manners, (i) The evolution of the effective action at its minimum, p = p gr , determines 
mainly the ground state properties, i.e. the quasiparticle wave functions ^ n and spectrum E n . We shall be looking 
for the generalization of the HF and the KS schemes in this way. (ii) The evolution of the functional dependence away 
from the minimum determines the different potentials parameterizing the increase of the effective action for p ^ p gr . 
We shall not assume that the 'excitation' p — p gT is small, 
(i) The minimum of the density functional is reached at 

n 1 ( n * n §C[p ^ fi+r x2e *\ in* a 5C[p ^ ] \ (m) 

^-T^g^fV - g -—) = { 1 + G -—)-{ p - g -—)- (85) 

Making use of the definition (f76|) and the relation 
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dx 



r[p] 



A 2 e 2 1 



= d\C[p gr ] - d\p* ■ 1 



|P=Pgr 



2 C 2 



A 2 e 



• G • 



2^2 



A 2 e 
A 



SC[p gI ] 



1 /<5C[ Pg 



2\ Sp r A 
with d\C[p gI ] = d\C[p]\ p=Psi , we find at the minimum 



G- 



A 



G-dxG) ■ 



( $C[pg 



A 



• P 



(86) 



^ oo 

Ae 2 Tr£ • - = -Ae 2 ^(-l)"Tr Q ■ (c^[ Pgr ] ■ Q 



=\ n 1 
A 



G 



_l(5_C[p M A_ . AV 
2 ^ <$p P A 

A 2 e 2 - A 2 e 2 
+9a/0* • —r— ■ G ■ — — ■ p* - d\p* 



G-dxG 



A 



f SC[ PgI ] 

SC[ Pg r] 

Sp 



p 



A 

dxC[p g r]. 



(87) 



This equation is in general inconsistent with the evolution equation imposed at p =/= p gI since the functions parame- 
terizing the functional, the propagators, are overdetermined. By following the strategy of turning the renormalization 
group equations into a partial resummation of the perturbation expansion we retain the modification of our effective 
action within the framework of the ansatz in the region where the approximation is supposed to be the best, in the 
ground state p — p gr . 

(ii) We shall use in this section the parameterization C = C q fl + C eX ct up to bi-local terms. The first term, chosen 
originally as 

- nsTrlogG- 1 + iTrlogP" 1 - ^Trlog-A 
and controlling mainly the zero point quantum fluctuations, will be replaced by the simplified form 



1, 



CqflM = -^TrlogG- 1 + ^Trlog 



-Tr D ■ a Tr a • D ■ a • D 

2 4 



(89) 



according to the assumption about bilocality. The photon propagator T> is constructed by allowing the photon 
self-energy to become a local function of the density, 



VZl = D al + ^(O, Pa)Sa,b = ["A + a{<X, Pa )]5 a , b + E^ 



(90) 



as assumed in the KS scheme ||. The explicit dependence of the self-energy on the position is due to the localized 
single-particle states in the ground state and can be omitted in the homogeneous electron gas. The second term of C 
is supposed to govern the excitations and will be the sum of an ultra-local and a bi-local contribution, 



f °° i 

CoxctW = / U(x, Px )+ V — — 



-,P -7 



(91) 



We carry out the replacement e — > Ae and introduce A-dependent single-particle wave functions W niX in these ex- 
pressions. For the sake of simplicity, we do not show explicitly the A-dependence generated in this manner. We 
find 



C^l=-\d p a{a,p a )D a , h d p a{b,p b )+ £ -^p n A 



n (n+l,m+l) m 
Pb 



+S a , b [d 2 p U(a,p a ) + lD atb d 2 p a(a,p a )-^d 2 p a(a,p a ) ( D aiC a(c, p c ) + £ -^p n a [ 7 & 

V - J c n,m=0 U - m - Jc 



+2,m) pin 



(92) 



in this parameterization, with D a b = D a ^D bta = (D a ^) 
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B. Evolution equation for the density Green function 



We consider first the evolution equation (|83|) for arbitrary densities in an /1C2 scheme for the effective action of 
Eq. (82). The potential C[p] — f U(a, p a ) is the local fx part of the approximation. Taking the second functional 
derivatives of both sides of Eq. ( p3|) and keeping the terms up to the bi-local ones (n = 2) on the r.h.s, one obtains: 



1 



1 



dx(Q-%,b = 2Ae 2 - 5 a . b - d x d 2 U(a, p a )8 a , b + Xe 2 (g ■ - ■ g) a , a d 4 p U{a, p a )S a , b 



-Xe 2 



-Xe 2 



g. d 2 pU .g._.g 



g.-.g.dp-g 



G--r-G) Q a ,b + [g--r-g) g b . a 



b , a 



a . b 



d p U{a,p a )S a ,b 
d 3 p U(a lPa )dp(b, Pb ). 



(93) 



Separating the terms a — b and a ^ b on both sides and exploiting the symmetry g ab — g b ^ a the following evolution 
equations are obtained for the particle-hole Green function and the local potential 



dx 



{Q- l )a t a-&Ma t p a ) 



2Xe 2 



1 



1 



+ \e"(g---g) a , a dtu(a, Pa ) 



la— b 



2Xe 2 [g-d 2 u-g---g) d 4 p u( a , Pa ) -2Xe 2 [g ■ - -g \ g a , a {9;u(a,p a )) 



<>.\(<j l ),,.,, = --2x< 2 [g-^-gj g a , b d 3 p u(a,p a )d 3 p u(b,p b ). 



(94) 



The numerical integration of these equations with the initial condition obtained by diagonalising the Hamiltonian 
for a non-interacting system represents an algorithm whose computational requirement depends on the space-time 
resolution of the Green function g required. The density p* can be determined by integrating the first functional 
derivative of Eq. evaluated at p = p* in A, 



1 



d xP * a = {d x d p u ■ Q) a - \e 2 I dp{b, Pb )(g---g\ h g b>a . 



(95) 



C. Evolution of the ground state (/a) 



The numerically more involved but more illuminating truncation scheme is where the quasi-particle wave functions 
are introduced in parametrising the effective action. Since the quasi-particles are well defined for excitations around 
the ground state we consider the evolution equation corresponding to these schemes at the ground state only. Since 



Tr 



g-^ + G-D 



00 , 

= -2TrG -r-^2[G- 



n=0 



, 2 2\ 2 ™+ 1 

X z e z N 



(96) 



and 



TrG" 1 ■ d\G — (3 f 

J X 



JV 



- E "-^ A *-» + V * Kt ) E ^ (*«^) |x=x' 



n=l 



n=l 



(97) 



the expressions , ([H]) and the relation ( |l5|) give the evolution equation as a generalized HF equation for orthogonal 
wave functions, 



N 

„ — 1 J X 



(98) 



|p=Pg 
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imposed together with the auxiliary condition 



dE„ 



= 



|P=Pgr 



where the generalized HF functional 

H[V*,%E,p] = n iiee [^*,^,E] +H Coulomb + H photon [^*,^,p\ + H int [V*,V,p], 

is given by 



(99) 



(100) 



N 



Coulomb fiT, * 



z 2 A 2 e 2 X 



2m 



n=l 

^photon^. = ^iHo g ^_A + o-(p)+A a e 2 G) , 

N 1 

l , n =i ,/x 'y 1 ■ y| 



[ ' ,PJ / Sp 2P\ 5 P 9 A 



:SC[p] A 2 e 2 , 



(101) 



and the inhomogeneous term is 

oo , 

pR,[V*,V, E] = 2Ae 2 TrG ■ — • ^ f G ■ 

71=0 ^ 

+Xe 2 TrD G ■ D ■ cr(p gr ) 
. A 2 e 2 



, 2 2 \ 2n+l 

A e N 



A 



OO r 1 

Ae 2 ^(~irTr Q ■ {c^[ Pgl \ • (?) 



1 - D ■ <j(p gr ) 



-TrD 



1 f SC[ Ps 

2 ^ Sp 



P ■ 



d 2Xe^ & ( SC[p SI } 
A I Sp 



1 - cr(pg r ) • £> 

A 2 e 2 



A 

d\cr{Pgr) 

d\C cxc t [Pgr] 



(102) 



with d\a[ PgI ] = d\<j(p)\ p = Psi . Note that the auxiliary condition ( j99| ) does not determine the value of E n . Instead, the 
functions g n ,n'(t) and the constants E n appearing in the particle-hole propagator are fixed by the dynamical relation 
( |l2| ) for a given set of wave functions \& n . It is worthwhile noting that similarly to the usual HF method the spectrum 
E n plays the role of Lagrange multipliers in the evolution equation in order to arrive at normalized wave functions. 

Another quantity not determined by Eq. ( ]9^ ) is C cxc t[p gr ]- Different choices for this A-dependent number lead to 
different wave functions ^ n and spectrum E n , different quasi-particles in short. The choice of the basis used in the 
expansion (Q) of the electron field is a central issue of any approximation scheme. The guiding rule to define the 
quasi-particles is to minimize their residual interactions, the higher order perturbative contributions. Since we are not 
explicitly computing the latter we turn to another, hopefully closely related principle for choosing the quasi-particles, 
the simplicity of the evolution equation, !Z[i>*, 9, E] = 0. The wave functions ^ n and the spectrum E n can then be 
determined by the HF strategy by minimizing Ti[^!* , ^ , E, p gr ] for normalized wave functions at each value of A, i.e. 
by imposing 



6H[**,*,E,p] 
5^* 



|P=Pgr 



6H[**,V,E, p] 



= 0. 



(103) 



|P=Pgr 



One can construct the renormalization group flow with the accuracy of 0(5\) when the wave functions at A + SX are 
found by minimizing the functional TL[^* , "J, E, p gr ] computed at A. The initial value for p gr is p gr (A = 0) = p , 
the ground-state density of the non-interacting system. 

The simplest truncation is to retain the terms C(p 2 ) and 0(e 2 ) of the density functional which gives C[p] =const. 
and cr(p) = according to the leading order expression (fTq), a C2 scheme. The dynamical problem, the evolution of 
the matrix elements ( |12|) does not appear at this level and the remaining static (time independent) evolution equation 
realizes the usual HF scheme following what happens when the interaction strength e in H2 is turned on gradually. 
The evolution offers in this manner the possibility of approaching the true HF minimum in infinitesimal, perturbative 
steps. Our original assumption about the smoothness of the e-dependence requires in the one-loop approximation the 
absence of local minima of the HF energy along the path < A < 1 which would prevent us to find the true minimum. 
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The terms 0(p i ) and C(e 4 ) ignored in the HF equation represent the multi-body forces and the contributions of the 
vacant HF orbits. But even the truncation C =const. and cr(p) — gives an already improved HF scheme where the 
exchange interaction coming from the second line of Eq. (101) and the energy of the zero point fluctuations when 
72. = is solve d co ntain a partial resummation of the perturbation expansion in e. The dynamical HF functional 
given by Eq. ( (lOOj ) implies the time-dependence of the propagator due to the intermediate states appearing in the 
higher orders of the perturbation expansion. 

We find it remarkable that the main feature of the HF scheme, the determination of the structure of the quasi- 
particles making up the ground state by minimizing a functional, can be maintained when higher order radiative 
corrections are retained. But there are complications which can easiest be understood by adopting the language of the 
perturbation expansion and recalling that the higher order contributions contain excitations around the unperturbed 
ground state. These contributions correspond to excitations and de-excitations of the ground state and imply time- 
dependence. As a result, despite that the original HF functional used in the determination of the ground state is static 
and refers to the spatial dependence of the wave functions only, the one given by Eq. ( |100| ) involves dynamical, time- 
dependent quantities. Furthermore, the construction of the intermediate excited states appearing in the higher orders 
of the perturbation expansion is achieved by the simultaneous solution of the evolution equation for the 'potential 
energy' of the density functional which is mainly responsible for the excitations. 



D. Evolution of the potential energy (/i) 

The generalized HF equation does not determine the evolution of the potentials U(p) and a(p) and that of the 
matrix ^(™> m ). In order to find them, we have to consider the effective action in its full bi- local form, 



~(p - P *) ■ dxg- 1 ■ ( P - P *) - d xP * ■ g- 1 ■ (p - P *) + I d x u(a, Pa ) 



00 1 ( 1 11 

V — —p n ■ dxl {n ' m) ■ p m + d x -n.TrlogG- 1 + -TrW-Z^A] + -Tr D ■ a Tr a ■ D ■ a ■ D 

— ' 2n\m\ » •>■-■> i 



= -Ae 2 | Tr-^ • Q - Tr-^ • Q ■ C {2) ■ Q + Tr-^ • Q ■ C (2) • Q ■ C (2) • Q 



(104) 



We introduce the parameterization p = p gr + r\ and use the identity 



\ 2 e 2 



\p=P S i+V 



2^2 



A 2 e 



— P-A-P 



+ dx(C[p gr + V }-C[ PgI }) 



|P=Pgr 
1 



+V ■ d\Q 1 • (p gr - p*) + -f] ■ d x G 1 • r) - d\p* - G 1 - r) 



(105) 



to write the difference of the evolution equations imposed at 77 ^ and rj = as 
Ae 2 f>irTr \g ■ [(c^[ PgI ] ■ of - (c^[p gI + rj\ ■ g)') ■ ±\ 

n— 1 V J 



d x (c[ Pgl . + v }- c[ PgI ]) + - v ■ dxG- 1 - v + v - dxG- 1 ■ (p gr - p*) - dxp* ■ Q- 1 ■ v- 



(106) 



Note that there is no term with dx Pe r in this expression because the substitution p — > p gI + rj is made after the 
evaluat ion o f the derivative 9a in Eq. (105). As mentioned earlier, the inequality \r)\ < p gr is not required. 

Eq. ( 106 ) together with the evolution of the ground state will be studied below in a more detailed manner in the 
ultra-local and the bi-local levels for the density fluctuation 77. We shall consider the potentials first and the results 
obtained will be used to find the explicit evolution of the ground state. 



E. Ultra-local approximation (/1C2) 

The scheme C2 is based on the ansatz (|2|) with C[p] = To being a Ae dependent but p-independent quantity. The 
effective actions of such form make up the space A C2 . This scheme will be improved by introducing an ultra- local 
part containing the functions U(x,p x ) and a(x,p x ) which have no linear piece in the Taylor expansion around p = 0. 
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This is achieved b y kee ping the ultra-local pieces without linear term only in Eqs. (|89j) and (|9l|). The ultra- local 
component of Eq. (|l06| ) restricted into the orthogonal subspace ^ fl Ajr 2 can be written as 



(107) 



for the function 



^a(A, p) = ^^a(a, p) + U (a, p), 



(108) 



where D a , a is rendered finite by the UV cut-off, 



/a(A)=Ae 2 [g 



1 

A 



(109) 



and the notation / gr (f?) = /(pgr + v) ~ /(Pgr) was introduced to denote functions of the density measured from the 
ground state value. 

Two remarks are in order concerning Eq. (107). The parameter A controls the amplitude of quantum fluctuations in 
the internal space and therefore it is rather natural that the evolution equation is a diff usion equation in the internal 
space with A as time. Furthermore, this equation is linear for the combination (108) and does not determine the 
potential U(p) and the function a(p) independently. 

By assuming the initial condition er = U = ip = we have no evolution because Eq. (107) is linear and homogeneous. 
The constrained bi-local form of the effective action is sufficiently general and further renormalization takes place only 
due to the three-body interactions, at the tri- local level. 

The evolution of the ground state energy is rather simple since = d 2 p ip = 0. In fact, Tq is found by integrating 
out the differential equation 



d x T = 2Ae 2 TrG 



1 

A 



OO 

E 

n=0 



G ■ 



A 2 e 2 
A 



2n+l 



1 



A 2 e 2 



A 



A 2 e 2 



A 



• P 



(110) 



wher e th e quasiparticle wave functions ^ on the right hand side are obtained by solving the generalized HF equation 
Eq. (102). The leading order contribution of the photon determinant gives the HF exchange term and its remaining 
contributions together with the terms of 0(e 4 ) in To represent the corrections to the traditional HF result. 

It is instructive to follow what happens if the evolution equation is allowed to adjust the linear part of the local 
potential what has been fixed in the computation above by the first term of the ansatz (^). One can see that the 
ultra-local part of the effective action is then modified in such a manner that the density in the ground state as given 
by Eq. (80) involves the partially resummed particle- hole propagator. The improvement is that we have a partial 
resummation of the leading order self-energy insertion. This is in agreement with the observation that the integration 
of the leading order renormalization group equation resums the leading order self-energy insertions in the propagators 
p8[ , J32| . The weakness of this solution compared with the HF level is that the non-interacting quasi-particles are 
used in the construction of the 'free' particle-hole propagator. This happens because we did not let the bi-local part 
of the evolution equation to correct for the correlations in the ground state and left the ultra-local part alone to 
determine the ultra- local potential and with it the location of the minimum in the space of density configurations p x . 



F. Bi-local approximation (/2) 



The bi-local part of the evolution equation (106) reads as 



E 



Ac 



d ^Ta 



(n,m) 



2n\m\ 



[(Pgr + V)a (Pgr + V)T ~ P&uP&bl ~ J d \ £>a,b {v(Pa + Va)v(pb + Vb) - <j(p a )o{pb)) 



P=P S r 



[ -r Gd,aGb,c\ y~] 



(n+l,m+l) 
Ta.b 



[(Pgr + V)2(Psr + V)b ~ PeraP&b 



-^D a , b [d v a(a, p gla + r} a )d v a{b, p gr6 + %) - d p cr(a, p gla )d p a{b, p glb )} 

-[d^1p a (X, p gro + 77 a )9 2 '0 & (A, p grb + T] b ) - dptp a (\, p gIa )dp^ b {\ Pgrb)]Ga,b 
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--^>a,b J D a , e [d%a(a, pgr a + r] a )a(e, p gre + r] e ) - dp(r(a, p gra )(r(e, p gre )] 

oo « (n+2.m) ~| -. 

+<y -f E / + ^"(fe + V)T - PgraPg^e] [ " ^A^W (HI) 

n,rn=0 ,/e J 

After taking into account the ultra-local solution = ?/> = and projecting out the bi-linear p-dependence present 
already in the 0(p 2 ) C2-part of the effective action this equation reduces to 

E -2M~Kp* + ^(Ps* + V)r - PlaP&b] 
( (n+l.m+l) 

= Ae 2 / - &, a ft, c J] ^ [( Pgr + 77)" (Pgr + »?)r - PgraPg^] 

Jc,d |^ nimi 

P (n+2,m) ~] 

+S a , b J2 / 2!i£ r^[(P gr + ^)"(Pgr + V)T ~ PlaP^re] (H2) 
Je Tl.Vf},. I 

where the summation over n, m = 2, • • • , oo, (n, m) = (1, 2) and (2, 1). Since this is a linear, homogeneous differential 
equation and the initial condition is 'y( n > Tn ) = 0, the solution is identically vanishing. Together with the similarly 
negative result of the ultra-local evolution one concludes that corrections to the scheme ci are generated at three- 
body forces only. 



G. Exchange energy in the j\Ci approximation 

It was noted above that the exchange energy comes from the photon fluctuation determinant. The evolution in 
the strength of the Coulomb interaction 'turns' the photon dynamics gradually on and should by itself generate the 
photon contributions to the effective action. Thereby, one wonders if the evolution can generate most of the exchange 
interaction without our choice of the ansatz for the density functional, i.e. by neglecting the particle-hole contribution 
to the photon propagator, D — > — 1/A. Let us simplify matters by restricting ourselves to a constrained bi-local 
approximation, f( m > n ') = 0. 

We note first that the basic idea of the K S sche me, the trade off the exchange energy for potential energy is nicely 
realized by the approximation /i of Section VIE. In fact, the evolution equation in this approximation is invariant 
with respect to the transformation 



U (a, p) -> U(a, p) + e(a, p), cr(a, p) -> a {a, p) 



D„ 



(a,p). 



(113) 



In the absence of the density dependence in the exchange interaction, as noted in section |VD| the success of the KS 
scheme can be ascribed to the efficiency of the ansatz where a local potential is kept, U(p) ^ in the absence of 
exchange interactions, D^ 1 = —A, 



The evolution equation, 



1 - A 2 p 2 1 

-(p-p*).g-i.( p -p*)-±— p .-. p . 



^TrlogG^+ro 



(114) 



~d xP * ■ g- 1 ■ (p - p*) + d x T - n^ATrlogG- 1 
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(115) 



considered at the minimum of the effective action, p — p gr simplifies to 



dxp* 



A 2 



■ p* - n s dxTi log G- 1 + dxT = -Ae 2 Tr Q • A + dxp* ■ ■ G ■ 
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(116) 
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because = 0. We can write this equation as 



N 

E 



"A^v FvtvI H **n,x- 



= 



together with the auxiliary condition (|99|). The generalized HF functional 

Wks[**, #, £?, P~] = W free [**, *, E] + H^ lomh [9*, % p] + H£i[V*, 9], 
is quadratic in the wave functions in the leading order of e, 

N 



*, Pi = n.AV / ^ *; x ^ n , x — J: 

■/x,y n=1 47T|X 

and To is given by the differential equation 

dxT 
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(118) 

(119) 
(120) 

(121) 



together with the initial condition To (A = 0) = 0. 

The solution of the evolution equation is given by the wave functions at the minimum of the functional HkSj the 
solution of the HF equations (103) together wit h th e auxiliary condition (|99|) obtained by the replacement 7Yhf —> Hks- 
The first term on the right hand side of Eq. (121) is canceled in Eq. ( |110[ ) against the leading order contribution of 
the photon determinant, the HF exchange term. In the present ansatz there is no exchange term but the evolution 
equation generates one since 



4^4 

2 A 



G(e 4 



In DFT one separates the kinetic energy T[p] and the direct Coulomb contribution by writing 



E v [p}=T[p] 



"47r|x — y 



■p y + E xc [p], 



(122) 



(123) 



and letting E xc [p] for the exchange contributions and radiative corrections. The effective action r free [/9] introduced 
above for the non-interacting fermion problem gives according to Eqs. (E3h and 



Tip] = lim ir free [pl = lim - 



^ ee + ^(p-p*)-G- x -(p-p*) 



+ o(( P -p*y) 



(124) 



The exchange and correlation functional r xc [p], the dynamical generalization of E xc [p], introduced in Eq. ( |123] ) for 
static density only, is defined as 



T xc [p]=T[p}-T^[p} + e Y p-^-p. 



In the approximation f\Ci we find 



r xc [p] = \{p - P*) ■ (g- 1 + £). (p - p*) - ^TrlogG-^W + r , 



(125) 



(126) 



c.f. Eq. (114), where G^ is the fermion propagator in the non- interacting system. 

It is important to note that although the effective action formalism is in contradiction with the basic assumption 
of the KS scheme this latter can be realized in a given truncation of the density functional. In fact, it has been 
noted that the exchange correlation, the one-loop contribution to the photon path integral, represented by To in Eq. 
(126) is independent of the density. But an artificially introduced density-dependence of the photon self-energy can 
be traded against a local density-dependent potential in the approximation f\ and opens the possibility of realizing 
the KS strategy. 
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VII. SUMMARY 



A renormalization group motivated method is introduced to perform the resummation of the perturbation expansion 
to the effective action for the density, the density functional, by solving a functional differential equation for the 
effective action. It is shown that this evolution equation which corresponds to the gradual increase of the strength of 
the Coulomb interaction generates the electrostatic energy of the charged particle system. In practical applications one 
has to project this equation on a restricted subspace of the effective action functionals and a systematic approximation 
scheme is suggested by truncating the functionals according to their multi-local structure. 

The simplest non-trivial truncation level reproduces the traditional HF scheme. The radiative corrections on the 
bi- local level where two-body correlations are retained only are captured by our generalized HF equations. The more 
involved approximation schemes provide a systematic generalization of the HF method by keeping the three-body and 
higher correlations. The most important complexity of these improved schemes comes from the fact that the higher 
order corrections of the perturbation expansion include contributions from virtual excitations which in turn renders 
the functionals in consideration dynamical, i.e. containing the time-dependence of the density as the independent 
variable. The generalized HF equation determines the wave functions of the quasi-particles making up the ground 
state by simultaneously minimizing the dynamical HF functional and solving the evolution equation for other functions 
characterizing the effective action. 

Generator functionals offer a formal way to separate the direct and the exchange energy contributions according 
to the level they enter in the multi-local expansion of the functionals or the saddle point expansion of the photon 
path integral. The direct and the exchange terms are responsible for the local, classical and the bi-local, one-loop 
contributions, respectively. 

The exchange term identified in this manner is density independent, in contrast to the basic assumption of the 
KS scheme. But it is shown that the claimed density dependence of the correlation energy can be traded into a 
local potential, independent of the structure of the quasi-particles in a given approximation of the density functional. 
The evolution does in fact generate in this case the correct exchange contribution to the ground state energy and 
the corresponding generalized HF equation gives an improved KS scheme. An explicit expression is obtained for the 
exchange-correlation functional. 

Finally, we mention two possibilities concerning the choice of the quantum state considered. According to the 
Wick's theorem, the expectation value (>I'|A| , I'} of the operator A depends on the state by the choice of the 
non-interacting path integral, the free propagator, and the Feynman rules which give account the interactions in a 
manner independent of \^>). In our scheme the interactions determine the right hand side of the evolution equation 
( |55| ) and the state |\&) is reflected in the parametrization of the electron propagator. One possibility is to generalize 
the ground state which according to the ansatz ( [l2|) corresponds to a single Slater determinant. By the introduction 
of an over-complete, non-orthogonal basis and the parameterization G XlX ' = 5Zn°n'=i 9n,n' (t ~ t')^n,x^n' x' we 
can mix several Slater determinants in the ground state. Another remark is about the issue of excited states. Suppose 
that the generalized HF equations (103) have been solved, a set of wave functions has been obtained and N of 
them, with the lowest 'eigenvalues' E n have been used to construct the ground state. The repetition of these steps 
when some of these wave functions are replaced with others which correspond to higher E n in the construction of the 
electron propagator (jL2|) leads to the systematic construction of excited states. 

We believe that the general setting of this computational algorithm for the density functional presented here, the 
gradual turning up the Coulomb interaction and a systematical approximation scheme by considering richer functional 
classes in solving the evolution equation is useful for strongly correlated electronic systems in Solid State or Co ndens ed 
Matter Physics, as well as in Quantum Chemistry. The numerical implementation of the scheme of Section VI B to 
test our ideas is in progress. 
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APPENDIX A: EFFECTIVE ACTION 



In this Appendix we briefly review the interpretation of the effective action for a generic scalar field theory. 
Suppose that the basic field variable of the system considered is (f> x and the connected Green functions are given 
by the generator functional W[j] expressed in terms of the path integral 

e wm = J v^-^W+L (Al) 

where Se[4>] denotes the imaginary time, Euclidean action. The effective action r[$] is defined by the Legendre 
transformation 
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m = -W[j] + [ $J X , <S> X = (A2) 

Jx °Jx 



(A3) 



The inverse variable transformation, 

shows that the minimum of the effective action gives the free energy in the absence of the source. 
In order to see what happens with non-vanishing source, we consider 

e" r [*] = max / VW^^L^^ (A4) 
ix J 

where the source j x is chosen by maximizing the right hand side. The field <j) x fluctuates in the path integral around 
the maximum of the integrand which satisfies the classical equation of motion, 

0(f>x 

By ignoring the fluctuations we have 



= mm I S E [4>] +J(4>x- <£x)jx) = min ^S E [<t> ci ] + J - <f> x )j x j , (A6) 

where we should bear in mind that </> cl depends on the source j. The minimization with respect j yields $ — cf> cl and 
r[$] = .Se^] showing that on the tree-level where the quantum/thermal fluctuations are ignored the effective action 
reproduces the classical one. 

In order to assess the role of fluctuations, let us assume that the system is placed into a large but finite spatial 
volume, introduce a basis 4>n,x, J 4>n,x4>m,x = fin.rm fi = 1,2,... for a given <& x where the element n = 1 agrees with 
&x, &x = Ccj)i tX and write the path integral as integration over the coefficient c„ in the expansion <p x = J2 n c « 



-nyn,x j 



r I*l = max f[ j d c „ e - SE[cl+(ci - C)jl+ £n 2 C "^, (A7) 



max 

3 k 



n=l 



We define the norm for the field configurations as \\<j>\\ 2 — J <t> x - Clearly |C| = indicating that the typical values 
in the path integral are c„ = 0(|| < f > ||°) for n — 2,3, .. . and C\ takes values around c\ — 0(||$||) with fluctuations 
5ci = 0(11$ ||°). We see in this manner that the relative fluctuations of the modes c„ in the path integral are weakly 
influenced by the source in directions orthogonal to the desired field expectation value $> x and are suppressed along 
the direction of Q x as 0(||<I>|| _1 ). As j is chosen such that & x extends over a larger volume the effective action tends 
to the free energy computed by integrating over configurations which have a fixed projection on <f>. 



APPENDIX B: FREE FERMIONS 



We re-derive in this Appendix the free energy of non-interacting fermions in the framework of Grassmannian path 
integration with special attention payed to the continuum limit, A< — » 0. 



1. Grassmannian path integral 

The fermionic amplitudes are usually obtained by integrating over Grassmann variables ]29| tpiip2 + ip2ipi = 0. The 
value of the integral J Uij)e~ s ^ is defined in such a manner that the equation of motion is satisfied on the average, 
/ V%lje- s ^6S[ip]/dijj = 0. This amounts to the rule 

/ *»m - ^ (bi) 

since any function f(ip) of a Grassmann variable ij) can be written as f(ip) = 77i'0 + 7 72 where 771 and 772 are Grassmann 
variables. One introduces pairs of complex conjugate variables ijr k , ipk and the integral 



25 



1 = II / d^ k d^le^ Ak -^ e = BctA (B2) 



can be easily checked by expanding the exponential function of the integrand and applying the rule (p3l|). 

Consider now a two-level system whose wave function can be written as ^f(ijj') = a + bip^ = \a, b), where a and b 
are complex numbers. The complex conjugate will be given by = a* + b*-ip = (a, b\. The states |1, 0} and |0, 1} 

are orthogonal and of unit norm by using the scalar product 

(*|*) = / d^d^e - ^***^)*^). (B3) 



We are looking for the matrix elements of the time evolution operator 

{a f ,b f \U(t)\a t ,b t } = (a^le-W^KM. (B4) 

For this end we divide the time interval [0, t] into N parts, t = NAt, introduce a pair of variables ijjt, ipk at each 
division point, and inspect the integral 

In = ft / #^*}(W)**(^)e^ I -^ W, -*- l) -* lfA ** J *'- l] 

fe=0 

JV „ JV 

= n / d4d^* f ^ N+1 )^i(4)l[e-^ ■ e^-^At)^.^ (B5) 

k=0 J j=0 

by expanding the exponential functions. The non-vanishing contributions proportional with must pick up from 
the first exponential function. This contribution is independent of E and will be responsible for the propagation of the 
ground state. In fact, always a pair ipj4>j what comes from the expansion of this exponential and (0, l|f/(t)|l,0) = 0. 
The contributions containing bi must find from the second exponential, responsible for the propagation of the 
excited state. We find in this manner 

™-j±. ( (i -r t)K ?) (B6) 

in the basis 

|0> = |1,0>=(5), |l) = |0,l)=fj). (B7) 

Notice that in order to have the right sign for the matrix element (1, 0|[/(f)|l, 0) we have to impose = ~i J N, 
ip^_i = — V>jy The matrix elements of the imaginary time evolution operator are obtained in a similar manner, 

h = ft / #|# fe */(^ + i)^(^)e E ^ h ^ t( ^^- l) - EA ^ t ^- 11 . (B8) 

fc=0 



2. Continuum limit 

The naive continuum limit, At — ► suggests 

= J Vl^}V[^}e* fy^ 9 *-^* (B9) 

and 

(* f \e-W*\Vi) = / V[i,]V[^]e-L^ +E ^\ (B10) 



2G 



The resulting continuum expressions are simpler th an th e discrete equa tion s but this simplification is rather misleading. 
In fact, let us evaluate the continuum expression ( BIO ) by means of (B2). The eigenmodes of the quadratic form in 
the exponent are 



and the corresponding eigenvalue is 

= no - E, 

where u) = (2n + \)~k j 3 due to the anti-periodic boundary condition in time. We have in this manner 

k—oo 



Det[-3t-E]= Yl (E - (2n + l)inT) = ]J ((2k+l)inT) JJ 



n— — oo 



n— — oo 



(2n + l)wrT 



By means of the relation pT[ 



n 



(2£+l)iTT 



— cosh ■ 



we find 



EB 1 ep 
Det[-ft -E] = C{T) cosh = -C{T)e~ 



1 



-B/3 



(BH) 
(B12) 

(B13) 

(B14) 
(B15) 



One might interpret the factor exp(EB/2) as the result of some kind of the zero point fluctuations and interpret the 
rest as the partition function for the two level system. But this is unacceptable because the first equation shows 
that the s uppre ssion is actually symmetric in E — > — E. We ended up wit h s uch an unacceptable result because the 
spectrum ( B12 ) is for continuous time, when At = 0. The expressions (B5), (B8) were obtained for At > 0, when we 
have 



■jAt\ 



-E = \ c u +0{At). 



(B16) 



What went wrong in the naive continuum limit? The high frequency contributions to the functional determinant 
which are essential in the case of a first order differential operator are missed in A c . The strong sensitivity on the 
details of regulating the action at the scale At can be seen by moving the variable ipj in the first term of the exponent 
in (|Bg| ) by At. It requires a change of sign of the kinetic energy and replacing the term handling the excitation by its 
complex conjugate, 



It = 



N 



(B17) 



Another modification appearing in this formula is that the wave functions are ^(ip) — a + bip and the excited state is 
represented by ^(ip) = ip- 

The problem arises because the trajec tori es ar e sup posed to have continuous first derivatives in classical calculus, 
behind the manipulations leading from ( p3§| ) to (BIO). This assumption is incorrect in quantum physics. It is well 
known that the Brownian motion trajectories are nowhere differentiable. Such a singular structure is natural in the 
path integral solution of the diffusion problem, 



P(X, T) = J\f~ 



lim 



N 

n 



dxje 



(B18) 



since N 2 (xj 



_i) 2 /T 2 sa 2DN/T — > oo for the typical trajectories p3|. The real time analogy of this phenomenon 



is that the trajectories in the path integral representation for the solution of the Schrodinger equation are nowhere 
differentiable. This problem becomes more serious when the ord er of t he differential equation, the power of the 
velocity in the exponent of the path integral is decreased, cf. (B8), ( B18| ). Another source of concern is the formal 
nature of the Grassmann variables. They po sses s no n umer ical values and obey no concept of smoothness. These are 
the actual sources of the difference between (B8) and (B17). 
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The operator G introduced in Eq. (Q) will be replaced by the regulated expression 



x.n.x'n' 



At 



(Sn 



) + s 



n,n' + l 



2m 



Ax - M + V x c 



Sx,: 



(B19) 



where the notation t = nAt is used. In this manner the first t erm in the right hand side of Eq. (|7(]) is obtained. The 
loop integrals containing the propagators obtained from (B19) have no more unexpected cut-off effects and reproduce 
the usual continuum expressions for time scales longer than At. As a result we shall deliberately use the continuum 
propagators given by Eq. ( |l0| ) in our computation. 



APPENDIX C: LEGENDRE TRANSFORMATION 

When the Legendre transform is sought it might be advantageous to convert the problem into those of solving a 
differential equation. This analogy which gives the solution of Eq. (^3|) in the simplest manner is discussed in this 
Appendix. 



1. Legendre transformation and differential equations 



Consider the non-linear differential equation f(y') — ty' + y = for the function y(t) and the initial condition 
y' (to) = vo where f'(vo) = to and f{x) is a convex function, f"(x) > 0. It is obvious that the pair of functions y(t) 
and f(x) are Legendre transforms of each others. As an example consider the choice f(x) = ax 2 /2 + bx + c, vq = 
which leads to ay' 2 /2 + (b — t)y' + y = c and y'(b) = and the solution y(t) = t 2 /2a - tb/a + b 2 /2a — c. 



In a similar manner, the Legendre transformation ( p3[ ) can be written as a functional partial differential equation 
and initial condition. 



W 



1 5T[p] 



i 5p 



n xr(l n 5r >°] 



Sp 



dp 



(CI) 



where <5W / [<7o]/ ( 5 cr = ipo- The quadratic generator functional W[a] — ^a-A-a + B- a + C gives the functional partial 
differential equation 



2 Sp 



Sp 



Sp 



Sp 



having the solution 



r M = ~\p ■ A' 1 ■ p - ip ■ A- 1 ■ B + ~B ■ A' 1 ■ B - C. 



(C2) 



(C3) 



2. Non-interacting electrons 

The functional Cauchy problem of the effective action of the non-interacting system is 



r froo [p] 



Trlog^-^^Up-^M, 



Sp ) 



Sp 



Sp 



0. 



p* 



Expanding the logarithm and using the Taylor expansion (|3q ) of p in powers of <r, one can rewrite the functional 
differential equation for r free [p] in the form: 



(C4) 



r frcc [p] = -^TrlogG- 1 - n s J2(l - -) f S Xl , X2 ,..., x „ 



8T boo [p} ST tlcc [p} ST bee [p] 



Spxi 5p X2 



$Px„ 



(C5) 



Looking for the solution in the form of the functional Taylor expansion (|43j), we can identify the expansion coefficients 
by taking the various functional derivatives of both sides of Eq. (C5) in increasing order. This yields just the result 
given in Eq. (E3). Thus, we recovered our previous result for the non-interacting electron gas. 
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3. Radiative corrections 



In the leading order one finds the functional (the right hand side of Eq. (71)) 



Wo + i (-J ■ A - iep* J • D ■ (-AJ - iep* 



1 e 2 1 - 

r J ■ A • J = Wo - —p* -D-p* -ip* -D-A-J+-J-G-J (C6) 
•'' 2 



for W[J] figuring in Eq. (CI) with J = which, according to (C2), reproduces j75|). 
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